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Dust yields in clumpy SN shells: SN 1987 A revisited 
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ABSTRACT 

We present a study of the effects of clumping on the emergent spectral energy distribution 
(SED) from dusty supernova (SN) shells illuminated by a diffuse radiation source distributed 
throughout the medium. These models are appropriate for Type II SNe older than a few hun- 
dred days, when the energy input is dominated by 7-rays from '"^^Co decay. 

The fully 3D radiation transport problem is solved using a Monte Carlo code, MOCASSIN, 
and we present a set of models aimed at investigating the sensitivity of the SEDs to various 
clumping parameters. We find that, contrary to the predictions of analytical prescriptions, the 
combination of an optical and IR observational data set is sufficient to constrain dust masses 
even in the case where optically thick clumps are present. 

Using both smoothly varying and clumped grain density distributions, we obtain new 
estimates for the mass of dust condens ed by the Type II SN 1987 A by fitting the optical 
and infrared spectrophotometric data of IWooden et al.l (Il993i) at two epochs (day 615 and 
day 775). When using amorphous carbon grains, our best fits to the observational data imply 
that about 2. 010^* of dust had condensed in the envelope of SN1987A by day 615 and 
between 2.0-10^'* and 4.2- 10^'^ by day 775. We find that the absence of a silicate emission 
or absorption feature in the observed mid-IR spectra implies that no more than 15% of the dust 
formed around SN 1987A can have been in the form of silicate particles. Our models require 
larger dust masses for the case of graphite grains, namely between 4.2- 10^"* and 6.6-10^^ M© 
at day 615 and between 4.5 •10^'' and 7.5 •10^'* M© at day 775. From our numerical models we 
derive dust masses for SN 1987 A that are comparable to previous analytic clumped graphite 
grain mass estimates, and at least two orders of magnitude below the 0.1-0.3 M© that have 
been predicted to condense as dust grains in primordial core collapse supernova ejecta. This 
low condensation efficiency for SN 1987 A is in contrast to the case of SN 2003gd, for which 
a dust condensation efficiency as large as 0. 12 has recently been estimated. 
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1 INTRODUCTION 



The production of dust in the ejecta of cor e-collapse super- 
novae (CC-SNe) is suppor t ed by observations (lLucvetal.1 [19891: 



Bouchet & Danziged [l993l: IWooden et all 1 19931 : lElmhamdi et al.l 
20031 : ISugerman et al while theoretical studies of dust con- 



densation in primordial CC-SNe have predicted th at 0.1-0.3 Mq 
of dust could be prod uced ( Kozasa, Hasegawa & Nomotd Il99ll : 
iTodini & Ferrarall200lh . Evidence for dust formation by (at least) 
some SNe comes from precise isotopic abundance ratio studies of 
grain inclusions found in meteorites, with many (including graphite 
and silicon carbide inclusions) exhibiting isotopic distributions that 



differ significantly from those found in the Sun and Earth (e.g. 
IClavton. Amari & Zinneil Il997l : iTravaglio et ai]|l999l) . However, 
direct evidence that SNe play a major role in the dust budget of 
galaxies is still sparse, since of the few supernovae that have been 
probed for dust production, most appear to have formed much less 
dust than predicted by models. For example, analyses of the pho- 
tometric and spectroscopic evolution of SNe 1987 A and 1999em 
from 400-800 days after out burst yield ed estima tes of a few x 10~^ 
M(7) of dust formed pe r SN jLucv et a l. 1989; W ooden et all 19931: 
lElmhamdi et al.l 12003), which translates into condensation effi- 
ciencies, defined here as (mass of refractory elements condensed 
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into dust)/(mass of refractory elements i n ejecta), of 10 ^ 
jLucv et alJl989l : IWooslev & Weaver! 19951) . In contrast, for SNe to 
have produced the majori ty of the >1 0^ Mp, of dust found in h igh 



igh- 
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redshift (z > 6) quasars iBertoldi et al.. 2 003 ; Maiolino et al. 200' 
iRobson et al1l2004l : I Hirashita et al.ll2005h. their condensation e ffi 
ciency must have been around ~0.2 jMorgan & Edmundsll2003h . 

In the cases of both SNe 1987 A and 1999em, the above in- 
vestigators all separately noted that the dust masses could be sig- 
nificantly larger if the dust was clumped rather than distributed ho- 
mogeneously throughout the ejecta. The presence of clumped dust 
is supported by the wav elength-independ e nt circumstellar extinc - 
tion seen in some SNe jLucv et ai] |l99lL lElmhamdi et d]|2003h . 
and by the theoretical expectation that Rayleigh-Taylor instabil- 
ities should form in post-shock ejecta fchevalier & Kleii] |l978l : 
[Herant & Woosley 1994). Sugerman et al. (2006) confirmed this 
hypothesis for dust formation within SN 2003gd, by demonstrating 
with three-dimensional Monte-Carlo radiative transfer (RT) mod- 
els that clumpy, inhomogeneous dust models can require up to an 
order of magnitude more mass, than implied by either analytic esti- 
mates or by homogeneous dust models, to produce a given optical 
extinction and mid-infrared (IR) excess. Sugerman et al. found that 
clumped models implied that up to 2 ■ lO"'^ A/0 of dust had formed 
in the ejecta of SN 2003gd by day 678, corresponding to a dust con- 
densation efficiency of up 0.12, which is interestingly close to the 
efficiency of ^^0.2 required to match high redshift dusty sources. 

Since the previous, rather low, dust mass estimates for 
SN 1987 A were based on analytical techniques, we were there- 
fore motivated to re-analyse the available observational data for 
this object using modelling techniques similar to those used by 
Sugerman et al. for SN 2003gd. In addition to allowing for 
clumped dust distributions, we also study the effects that the ge- 
ometry of the medium and the spatial distribution of the lumi- 
nosity source have on the emerging spectral energy distributions 
(SED) and, as a consequence, on dust masses determined by fit- 
ting RT models to observations. Radiative transfer in an inhomo- 
geneous medium has been the subject of a number of studies, in- 
cluding the simplest approach of a two-phase medium consisting 
of dense clumps embedded in a less dense interclump medium 
(ICM). These studies employed a combination of analytical tech- 
niques and Monte Carlo modelling with dust models of varying 
complexity, in order to study dust in different phases of the in- 
terstellar med ium (ISM) and in pro t ogalactic and galactic envi- 
ronm ents (e.g.'Natta & Panagia"l984':'Boisse|ll99dl:lNeu feldll99ll: 
Hobson & Scheuei 1993; Hobson & Padman 1993; Witt & GordonI 
1996; Gordon, Calz etti & Witllll997l ; IWolf. Fischer. & Pfaulll998l ; 
'Varosi & Dwek \99^ 



All the models presented in this work were performed using 
the 3D Monte Carlo RT code mocassin jErcolano et~d]|2003al ; 
lErcoIano. Barlow & Storevl2005i) , which can account for both (pri- 
mary and secondary) components of the radiation field in a fully 
self-consistent manner through absorption, re-emission and scatter- 
ing of photons. This is necessary if one wants to deal with the trans- 
fer through inhomogeneous clumpy distributions, both for the case 
of a diffuse illuminating source distributed through the medium or 
for the case of a central illuminating source. Section 2 describes the 
basic assumptions used for the modelling, while the results of our 
parameter investigation are presented and discussed in Section 3. 
In Section 4, we readdress the issue of the mass of dust formed 
by SN 1987 A, constraining our models primarily with the obser- 
vation s between ~0.3 /im and 30 /im published by IWooden et al.l 
( I199I hereafter W93). 




Figure 1. A 3D representation of the dust number density in clumpy mod- 
els. 



2 MODELLING STRATEGY 

Using both smoothly varying and clumped grain density distribu- 
tions, the 3D Monte Carlo r adiative transfer (RT) code MOCASSIN 
tercolano et al.l2003a|; Ercola no, Ba rlow & Storey 2005.) was used 
bv ISugerman et al.l ( l2006h "to derive the mass of dust condensed 
by the Type II SN 2003gd, by fitting optical (B, V, R, I) and 
Spitzer Space Telescope 3.6pm, 4.5pm, 5.8/im, 8.0/im, and 24pm 
fluxes observed at two epochs (days 499 and 678 after outburst). 
Their results clearly showed that analytical analyses of optical and 
IR data, as well as RT models with smoothly varying dust den- 
sity distributions, can underestimate the dust mass by an order 
of magnitude or more compared to models that allow for clump- 
ing. ?? provided evidence that optically thick clumps may exist 
in the envelope of SN 1987A, implying that the dust mass may 
be underestimated by SED-based smooth density distribution mod- 
els. These authors argued that amorphous part dust was required 
by the wavelength-dependent blueward shift of various emission 
lines, while clumps are required to explain the wavelength indepen- 
dent reddening effect of the dust. Successively, W93 summarised 
the evidence for significant clumping in the ejecta of SN 1987 A, 
which was provided by observations ranging from 7-ray and X-ray 
to infrared wavelengths. We have therefore revisited the observa- 
tional data for SN 1987 A, employing a similar modelling strat- 
egy to that used for SN 2003gd, i.e. we assume that dust con- 
densed in clumps within the ejecta, with the heating source dis- 
tributed within the remnant shell. Our models are appropriate for 
Type II SNe older than a few hundred days, when the energy in- 
put is dominated by 7-ray s from ^''Co decay jBouchet & Danzigeil 
I1993I W93). We consider two clumping models, designed to rep- 
resent the extreme cases of likely clumping behaviour; (I) the dif- 
fuse heating source is mixed with the dust grains, which are present 
in both the clumps and in the interclump medium (ICM); and (II) 
the diffuse illumination source is confined to the virtually dust- 
free ICM, with all dust grains confined to the clumps. The latter 
model appears the most physically plausible, given that the outer 
C, O, and Si-rich shells, where carbonaceous and silicate grains 
are expected to form, should never completely mix with the inner 
zone jArnett etaljl989l) . Only macroscopic mixing has 
been foun d in the clum py ejecta of the Type II supernova remnant 
Cas A (Do uvion. Lag age & Cesarsky 1999), suggesting that the el- 
ement layers were not homogeneously mixed. For completeness, 
we also consider a "smoot h", clump-free model. 

lMoselevetal.l ( ll989l) reported early IR observations (day 265 
and 267) and found no evidence for any emission from dust that 
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Figure 2. Amorphous carbon model fits to the day 615 spectrophotometry 
published by W93 (black line). Two smooth density models of day 615 with 
different bolometric luminosity are compared. The red-dotted line is for a 
model with L = 5.7 X IO^Lq model and the green-dashed line is for a 
model with L = 6.7 X IO^Lq. 



might have formed in the SN ejecta or from pre-existing dust 
in the surrounding medium. A variety of observations indicated 
that the onset of dust fo rmation in the e jecta of SN 1987A oc- 
curred at about day 530 jLucv et alj|l991i . W93). Given that the 
featureless lO-fim spectrum of SN 1987A after day 530 (W93; 
iRoche. Aitken & Smith|[l993h precludes a significant silicate con- 
tribution, the dust was assumed to be either amorphous carbon 
(optical co nstants from iHanned Il988h or graphite (optical con- 
stants from lDraine & Le3 ll984'). In view of theoretical predictions 
that supernova gr ain s ize distributions should be relatively lacking 
in large particles ^To dini & Ferrara 200 ih . we adopted a standard 
jMathis, Rumpl & Nordsiecktl 977. hereafter MRN) ISM grain size 
distribution, N(a) oc a~^'^, but with the grain radii a truncated at 
the upper end, such that amin = 0.005 nm and Omax = 0.05 fim. We 
found however that for the models presented in Section 3, the SEDs 
predicted by considering a full MRN distribution versus a truncated 
size distribution are quite similar and that the choice of either did 
not affect the general results of this paper. Following our parameter 
investigation, we determine best fits to the W93 day 615 and day 
775 spectra for all three dust configurations (smooth. Clumpy I and 
Clumpy II) and in Section 4 present revised dust mass estimates for 
SN I987A on days 615 and 775. 



2.1 The 3D mocassin code 

The 3D Monte Carlo photoionisation an d dust RT code MO- 
CASSIN was described and bench marked in jErcolano et al.ll2003al : 
lErcolano. Barlow & Storevll200 5V We used version 2.02.21 of the 
code in its dust-only RT mode, which allows a multi-grid resolution 
approach in order to increase the spatial resolution in the clumps. 
The radiati on field is described by discrete monochromatic packets 
of energy jAbbott & Lucvl[T985h . whose trajectories through the 
shell (mapped on multiple Cartesian grids) are characterised by ab- 
sorption, re-emission and scattering events, according to the local 
medium opacities and emissivities. To ensure that energy conser- 
vation is strictly enforced at each location, all energy packets are 
followed until they escape the shell, contributing to the emergent 
SED. 



The initial frequencies of the energy packets, emitted in this 
case by an extended diffuse source, are determined by sampling the 
probability density functions (PDF's) obtained from the blackbody 
distribution for the temperature of the source. The total bolomet- 
ric luminosity, L, integrated over the diffuse source equals that of 
an equivalent centrally-located point source and the energy carried 
by each individual energy packet in unit time is simply L/A^'pack, 
where A'^pack is the total number of energy packets employed in the 
simulation. 



2.2 The dust density distribution 

Our models consist of a spherical shell, with inner radius Rin, 
outer radius 7?out = Y-Rin and a p oc r^" density profile. 
Two basic dust density distributions are considered: one where the 
density varies smoothly as and the other with dense homoge- 
neous clumps embedded in a less dense ICM, with a density pro- 
file varying as r~^. The clumps have radius 5 ■ Rout and a vol- 
ume filling factor, /. The density contrast between the clumps and 
the ICM is defined by a = Pc(-Rin)/picm(-Rin), where pc(-Riii) 
and picm{Rin) are the clump and ICM dust densities at the in- 
ner edge of the shell. The absolute dust density in the clumps is 
therefore constant for all clumps, regardless of their radial dis- 
tance, meaning that the density contrast with the ICM is larger at 
larger shell radii. Our choice to employ this model is based upon 
the assumption that as early as a few hours a fter outburst, post- 
shock ejecta become Rayleigh-Taylor unstable l lChevalier & Kleij 
1978; Herant & Wooslev 1994), forming a clumpy distribution. 
The clumps would then simply move outward with the original 
imprint of density enhancement remaining intact. The clump po- 
sitions are assigned stochastically, with the probability of a given 
grid cell being occupied by a clump being proportional to r^'^. The 
total number of clumps in a given model is controlled by the vol- 
ume covering factor, /, and the clump size parameter, which is kept 
constant at a value of S = f?out/30. 

The smooth models use 61'^ cells to describe the SN shell, al- 
though we make use of the assumed symmetrical nature of the ob- 
ject and restrict our simulations to only one eighth of the total vol- 
ume, r eflecting the packets at the Cartesian planes jErcolano et al.l 
l2003bh . Figure [T] shows a 3D representation of the density distri- 
bution in the clumpy models, consisting of cross-sections taken 
along the positive domains of the orthogonal Cartesian planes. The 
clumpy models also use a mother grid of 61'^ cells to describe the 
ICM, whilst each individual clump is described by a subgrid of 5'^ 
cells which occupies the volume of a mother grid cell. We exper- 
imented with higher resolution grids and subgrids and found that 
they had negligible consequences on our final results. 

For clumpy model configurations I and II we first vary the 
dominant clump parameters and show the effects of these variations 
on the emergent SEDs and equivalent optical depths, defined in 
Section |23] below. 



2.3 The diffuse illuminating source and equivalent optical 
depths 

The illuminating source is adopted to have a total bolometric lumi- 
nosity L = 5.7 X lO^L© at day 615 (from Figure 3 of W93) radi- 
ating as a blackbody of TcS = 7000 K. We experimented with lower 
blackbody effective temperatures and came to the conclusion that 
the observations could be best fitted by a value of Teg = 7000 K, 
which may be interpreted as the superposition of the 5000-5500 K 
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Figure 3. Clumpy I models: dependence of SEDs (left panel) and equivalent optical depths (right panel) on the volume filling factor, /, of the clumps. In 
this and subsequent plots, SEDs (left) are plotted as dimensionless, distance- and luminosity-independent spectral shapes, A F\l J FxdX versus A in microns; 
equivalent optical depths (light) are defined by Equation (1), and may thus be negative if the escaping flux is greater than the incident energy, as is the case 
in the mid-IR. All models have the same = 5.6-lO~*M0 and a = 100. Colour coding: black (dotted) for / = 0.01, red (dashed) for / = 0.02, green 
(dash-dot) for / = 0.05 and blue (dash-dot-dot) for /=0.10. 
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Figure 4. Clumpy I models: dependence of SEDs (left panel) and equivalent optical depths (right panel) on the clump to ICM density enhancement factor, a. 
All models have the same Md = S.e-lO-^M© and / = 0.05. Colour coding: black (dotted) for a = \, red (dashed) for a = 10, green (dash-dot) for a = 100 
and blue (dash-dot-dot) for a=1000. 



optically thick 'hot' component and the H I bound-free continuum 
component identified in the W93 fits to the SN 1987A observa- 
tions (thei£^i£ure_^)^_We ahon t hat for the same epoch (day 
616) ISuntzeff & Bouchetl ( Il99(]|) and lBouchet & Danziged ( Il991l) 
found L = 6.6 X 10^ L© and L = 6.7 x 10^ Lq, respec- 
tively. However in this work, we adopt the lower luminosity value 
which is consistent with and returns the better fit to W93's data, 
as shown in Figure |2l where we compare the SEDs obtained by 
identical smooth density models of day 615, but with different 
bolometric luminosity. The red-dotted line is for a model with 
L — 5.7 X model and the green-dashed line is for a model 

with L = 6.7 X lO^Z/0. The input parameters for all models will 
be discussed in more details in the following sections. The diffuse 
source is distributed either throughout the medium, including in the 
clumps (the Clumpy I models, in which the diffuse heating source 
intensity follows the dust density); or located only in the ICM (the 



Clumpy II models), such that each radiating cell i emits a luminos- 
ity Li proportional to the local interclump medium density. 

In the case of a SN shell illuminated by a diffuse source of ra- 
diation distributed within the shell, optical depths calculated along 
a line of sight from the centre of the sphere have little observa- 
tional relevance. To account for the fact that the photons emitted 
from different locations within the shell will encounter very differ- 
ent opacities, we calculate the effective optical depth, rcfr(A), from 
the ratio of the emerging unattenuated flux at A for the case where 
there is no dust present, F^, to the emergent flux after attenuation 
by dust in the shell, F\: 

pO 

roff(A) = ln-f (I) 

This method is completely independent of geometry and therefore 
does not require averaging along a discrete number of lines of sight; 
furthermore it automatically takes into account the effects of pho- 
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Figure 5. Clumpy I models: dependence of SEDs (left panel) and equivalent optical depths (right panel) on the total dust mass Aid- All models have the same 
a = 100 and / = 0.05. Colour coding: black (dotted) for Md = 5.6- 10"'^ A/©, red (dashed) for Md = IMO-'^M©. 
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Figure 6. Clumpy II models: dependence of SEDs (left panel) and equivalent optical depths (right panel) on the volume filling factor, /, of the clumps. All 
models have the same A/d = 5.6-lO-"'Af0 and a = 100. Colour coding: black (dotted) for / = 0.01, red (dashed) for / = 0.02, green (dash-dot) for / = 0.05 
and blue (dash-dot-dot) for /=0.10. 



tons scattering off the grains, which increases their path lengths 
through the shell, thus increasing the probability of their being ab- 
sorbed by dust before escaping. 



3 RESULTS: THE EFFECT OF CLUMP PARAMETERS 
ON THE EMERGENT SEDS 

3.1 Clumpy I models: diffuse heating source and dust located 
in both the clumps and in the ICM 

In this section we present results obtained from clumpy models 
where both clumps and ICM radiate with luminosity Li propor- 
tional to the local medium density. The effects of varying the vol- 
ume filling factor of the clumps, /, the clump-to-ICM density ratio, 
Q, and the dust mass, Md on the emerging SEDs and effective opti- 
cal depths of the models are illustrated in figures [3]to|5] The SEDs 
are plotted as dimensionless, distance- and luminosity-independent 
spectral shapes, or relative flux, A F\/J F\ dX, versus A in microns. . 
Figure [3] shows that increasing the volume filling factor / of 



the clumps has dramatic effects on the effective optical depth and 
therefore on the emergent SEDs, for models with a = 100 and the 
same total dust mass A/d = 5.6- 10~''A/q, with the individual dust 
clumps each having a mass of lO"''Af0. When / = 0.05, each of 
these clumps has a centre-to-edge radial optical depth in the R- 
band of 1.22 for the amorp hous carbon models . The same clumps 
containing graphite grains ( iDraine&Le3ll984t) . rather than amor- 
phous carbon grains, would have a centre-to-edge radial optical 
depth in the R-band of 0.46. 

It is clear that an increase of the clump volume filling factor 
causes a decrease in the effective optical depth of the models, with 
less optical radiation being reprocessed at IR wavelengths. This 
somewhat counter-intuitive effect is in fact a direct result of the 
luminosity source being embedded in the clumps, so for constant 
Afd and a, an increase in / implies a decrease in the clump dust 
density and thus a reduction of the optical depth through each knot, 
facilitating the escape of photons emitted by a given knot. This is 
opposite to the behaviour of models which have a central source of 
radiation, or ones in which only the ICM radiates, as described in 
the next section. 
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Figure 7. Clumpy II models: dependence of SEDs (left panel) and equivalent optical depths (right panel) on the total dust mass A/d. All models have the same 
a = 100 and / = 0.05. Colour coding: red (dashed) for A/d,o = S.S-lO'^Af©, black (dotted) for = 0.1xA/d,o, green (dash-dot) for M^i = 10xA/d,o, 
blue (dash-dot-dot) for A/j = 50x Afd,o and cyan (long dash) for Afj = lOOx Af^.o- 



In figure|4]we show the SEDs (left panels) and the equivalent 
optical depths (right panels) for models with a constant / = 0.05 
and a total dust mass in all clumps of Aid = 5.6-lO~'*M0, for dif- 
ferent values of a. Variations in a do not appear to have large ef- 
fects on the emerging SEDs, except perhaps for the a = I model 
(red-dotted line), where there is no dust density enhancement in the 
clumps with respect to the ICM. We note, however, that the ICM 
dust density follows an r^^ dependency, whereas the density of 
the clumps stays constant with distance from the centre, implying 
that the clump to ICM density ratio increases as the square of the 
distance from the centre. 

The main point is that for a model where the intensity of the 
diffuse heating source follows the dust density, dust grains cannot 
be hidden in optically thick clumps. An increase in dust mass, e.g. 
a factor of two, will change the optical and IR SED dramatically, as 
shown in figure [5] where in the left- and right-hand panels we plot 
the emergent SEDs and equivalent optical depths for models where 
/ and a stay constant at values of 0.1 and 100, respectively, and 
the total dust mass Afd varies from ~5-10~''A/q to -MO-H/q. 
Finally, it is clear that for these models a degeneracy exists between 
/ and total dust mass A/d, which cannot be resolved solely by SED 
observations. 



3.2 Clumpy ii: radiation from dust-free ICM and dust 
confined to clumps only 

We now present the results obtained from clumpy models II having 
exactly the same clump distribution as for I, but this time assum- 
ing that dust is only present in the clumps and that the luminosity 
source is only embedded in the ICM. The parameter a tends to in- 
finity for these models, so we only investigate the effects of varying 
/ and A/d on the emerging SEDs and equivalent optical depths. 

The left and right panels of figure[6]show the SEDs and equiv- 
alent optical depths corresponding to models with constant total 
dust mass AU = S.S-lO-^Afe and / = 0.01 (black, dotted), 0.02 
(red, dashed), 0.05 (green, dash-dot) and 0.1 (blue, dash-dot-dot). 
The left and right panels of figure [7] show the SEDs and equiva- 
lent optical depths corresponding to models with constant / = 0. 1 
and Md between (0.35 - 3500) x lO'^Af©. For the reference 
model, with / = 0.1 and and Afd = 3.5-lO~'*Af0, each clump has 



a mass of 3.5- 10^ Af© and a centre-to-edge radial optical depth in 
the V-band of 1.33 and 0.41, for the amorphous carbon and graphite 
models, respectively. 

It is clear from figure |6] that / directly controls the effective 
optical extinction, i.e. the amount of radiation intercepted by the 
clumps and reprocessed at IR wavelengths. The behaviour of the 
SEDs and equivalent optical depths here is the same as for a cen- 
trally located source and opposite to that shown by the Clumpy I 
models where the heating radiation is emitted within both clumps 
and ICM. 

In the current case of a dust-free ICM with the dust confined 
to clumps, for optically thick clumps the equivalent optical depth 
in the visual region is governed only by the covering factor of 
the clumps, which determines the amount of radiation intercepted 
by the clumps. This would imply that the actual dust mass in the 
clumps would be unconstrained by the measured extinction in the 
optical and that virtually unlimited amounts of dust could be hidden 
in optically thick clumps. This is apparent in figure|7] where it can 
be seen that the models with 10, 50 and lOOx 3.510"''A/q show 
the same SED between ~0.3/im and ~2^m. However the quantity 
of dust in the clumps still has a significant effect on the grain tem- 
perature distribution and therefore on the appearance of the SED at 
IR wavelengths. The larger the dust mass for a given /, the cooler 
the average grain temperature and therefore the greater the shift of 
the reprocessed radiation peak to longer wavelengths. 

So, for the Clumpy II case the total dust mass is not uncon- 
strained, even when the clumps are optically thick. If both opti- 
cal and IR observational data are available, it appears that the best 
modelling strategy for optically thick clumps, when segregation be- 
tween the dust and the heating source is assumed, is (i) to start with 
a large dust mass (enough for the clumps to be optically thick in 
the visual region) and to vary / until a fit is achieved for the optical 
part of the SED; (ii) keeping / fixed, to vary the dust mass until a 
fit to the IR part of the SED is also achieved. If the dust mass is 
reduced to the point where the optical SED starts changing again, 
this implies that the clumps are becoming optically thin. This may 
be acceptable, or one may try a lower covering factor and restart 
the fitting process. 
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Table 1. Input parameters for SN 1987A dust shell models. The parameters specified in the top section of the table are common to all models. These include the 
total bolometric luminosity of the diffuse illumination source, L,, assumed to radiate as a blackbody of temperature T{BB), the inner shell radius, i?,;,!, and 
the index p of the power law describing the radial dependence of the density distribution. The carbon dust is either 100% amorphous carbon (optical constants 
from Manner 1988) or 100% graphite (optical constants from Draine & Lee 1984), with a standard MRN grain size distiibution truncated at the upper end 
(o-min = 0.005 fim; amax = 0.05 /im). Both smooth and clumpy models with a diffuse illumination source are considered. For the smooth models the dust 
number density decreases as . The Clumpy I models have radius-dependent clump-interclump dust number density enhancement factors, a, corresponding 
to the enhancement at the inner radius of the shell. Clumps are randomly distributed within the shell, according to a probability density function proportional 
to the square of the inverse of the radial distance. The clump volume filling factors, /, are also given for each model, as well as the total dust masses, M^, for 
both smooth and clumpy models with amoiphous carbon or graphite grain and the outer-to-inner shell radius ratio, Y = -Rout /Rin - The R-band extinctions. 
An, are also listed for each model, as well as the radial centre-to-edge extinction within an individual clump, Aji{c\ump). 
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4 SN 1987A REVISITED 

Taking into account the results presented above, we now derive re- 
vised dust mass estimates for SN 1987A. Our analysis is based on 
model fits to the 0.3-30 ^m spectrophotometry from days 615 and 
775 published in figures 1 and 2 of W93, which were digitised using 
the Dexter t ool pro vided by the NASA ADS, then dereddened using 
the lMathisl ( fT990h extinctio n law, scaled to fit th e extinction mea- 
sured toward SN 1987A in IScuderietal.Nl996 ^. Our models are 
also required to match the optical extinct ion Ar of the SN , which 
was first shown to in crease with time by Lucv et al] ( Il99 ll) . How- 
ev er, as discussed inlSugerman et all ( |2006|) . the method adopted 
bv lLucv et al.l ( ll99ll) to measure the dust extinction from the broad- 
band light curves is likely to have significantly underestimated the 
V and i?-band values. 

Our smooth density distribution and our Clumpy II models for 
day 615 (red dotted and blue dash dot lines in Figure [8] for amor- 
phous carbon and Figure [9] for graphite) both fit the IR and opti- 
cal SED of SN 1987A quite well - the smooth density distribution 
graphite model (Figure 8) provides a somewhat better fit to the 20- 
/xm spectrophotometry than the clumped models, but at the expense 
of under-predicting the flux level in the short wavelength optical re- 
gion. The effective extinction predicted by our Clumpy II models 
for day 615 is A(U) = 1.69, A(B) = 1.66, A(V) = 1.61, A(R) = 1.54 
and A(I) = 1.48, for amorphous carbon grains and A(U) = 1.64, 
A(B) = 1.54, A(V) = 1.44, A(R) = 1.38 and A(I) = 1.31, for graphite 
grains. Although our absolute values are larger than those reported 



by ILucv et ^ il99lh - see discussion in ISugerman et alj ^20061) 
ab out the extinction being underestimated by the method adopted 
by ILucv et il99lh -, the slo pe of the comput ed extinctions is 
in agreement with the results of lLucv et alj jl99lh . who noted that 
the empirical extinctions were rather flat in the optical. Our day 
615 Clumpy I models, whilst providing satisfactory fits to the IR 
SED, provides a poorer match in the optical region. For the day 
775 observations (Figure [TO] for amorphous carbon and Figure [TT] 
for graphite), the smooth and Clumpy I models produce too much 
extinction at the shorter optical wavelengths, particularly for the 
graphite case. 

All inputs for our models are listed in Table 1, where we also 
summarise our dust mass estimates, R-band model extinctions and 
the R-band radial centre-to-edge extinction within an individual 
clump. 

The day 615 smooth density distribution amorphous car- 
bon model requires a total dust mass of 2.O-1O-*M0, while 
the amorphous carbon dust masses derived from our Clumpy I 
and Clumpy II models are the same within the uncertainties, 
2.O-1O-*M0 and 2.2-lO-*M0, respectively. Our graphite grain 
models require approximately 3 times more mass to fit the SED 
of SN 1987A than do our amorphous carbon models, due to the 
significantly lower o ptical to mid-IR absorption coefficients of 
iDraine & Le3 ( 1 19841) graphite compared to Hanner (1988) amor- 
phous carbon (see Figure 112b . Our smooth density distribution 
graphite model requires a dust mass of 4.2-10-*Mq, a factor of 
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Figure 8. Amorphous carbon model fits to the day 615 spectrophotometry 
published by W93 (black line). All the models used an MRN size distri- 
bution truncated to maximum grain radii of 0.05 ^m. The red dotted line 
corresponds to a smooth density distribution model, the dashed green line 
to our Clumpy I model and the blue dot-dash line to our Clumpy II models. 
Model parameters are summarised in Table[T] 



Figure 10. Amorphous carbon model fits to the day 775 0.3-30-/im spec- 
trophotometry publishe d by W93 (b l ack lin e) and the day 793 50- and 
100-/^m photometry of i Harvev et al.l )l98^ (tiiangles). All models used 
an MRN size distribution truncated to a maximum grain radius of 0.05 /im. 
The red dotted line coiresponds to a smooth density distribution model, the 
dashed green line to our Clumpy I model and the blue dot-dash line to our 
Clumpy II models. Model parameters are summarised in Table[T] 
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Figure 9. Same as Fig.[8]but for graphite model fits. 




Wavelength [/im] 

Figure 11. Same as Fig.[lO]but for graphite models. 



10 higher than estimated by W93 for day 615 using their simple 
large-grain analytic model. Our Clumpy I and II graphite models 
both require ~6.5-lO~*M0, factor of two larger than derived by 
W93 for day 615 using an analytical clumpy model with graphite 
absorption coefficients. 

For day 775, our smooth density distribution amorphous 
carbon model requires a dust mass of 3.0-10-''Mq, while our 
Clumpy I and Clumpy II models require 2.0- 10~ M0 and 
4.2-lO~*M0, respectively. Our smooth, Clumpy I and Clumpy II 
graphite models (Fig. lilt require masses of 4.5-lO~*M0, 
6.5-lO"*M0 and 7.5-1O"''M0 respectively to fit the day 775 
SN 1987A data. The day 775 smooth density distribution amor- 
phous carbon and graphite models both provide better fits to the 
observed SEDs landwards of 20 /im than do the clumpy models, 
but both of the smooth density distribution models fail to match the 
observed flux level in the short wavelength optical region. The IR 
spectrum seems to be too sharply peaked at short wavelengths to 
be fitted with simple models. The residuals are also broader than 



typical atomic line widths, which could be interpreted as evidence 
for solid-state features in the dust -possibly stochastically heated- 
as discussed Bouchet et al., 2004. Our smooth density distribution 
graphite mass estimate is 7 times larger than W93's analytic large- 
grain estimate for this epoch, while our Clumpy I and Clumpy II 
graphite grain mass estimates are factors of 1.3-1.5 larger than 
W93's analytic clumpy model estimate of 5.0- lO^^M© for graphite 
grains. Dwek et al. (1992) also estimated a dust mass of 5-lO"*M0 
for the ejecta of SN 1987 A, based on fitting a 150 K modified 
blackbody to day 1 144 IR photometric data. Although their esti- 
mate is somewhat uncertain, due to the small number of photomet- 
ric points in the mid-IR region, it does indicate that the mass of 
dust in the ejecta of S N 1987A did no t increase significantly after 
day 775. Furthermore. iBouche 3 l2004h presented mid-IR imaging 
of SN 1987 A and reported, in particular, the detection of the SN 
ejecta in the N-band. They concluded that their measurements im- 
ply a dust temperature in the 90-100 K range, and a mass range 



Dust formation in clumpy SN shells 9 



MRN; a(n- 



nax) = 0. 005:0. 05/irT 




Jftot) qraph {Draine & Lee, 1S_ , 
Qtabs) graph (Draine & Lee, 1984 

Q(tot) omC (Hanner,19^_ 

" Q(abs) omC (Hanner,1988) 

Q(tot) est Sil (Draine & Lee, 1984) 

Q(obs) osl Sil (Droine & Lee, 1984) 

1 0"'° I , \ , \ , \ 




Figure 12. Absorption and tota l extinction coeffi cients (Q) of lPraine & Led 
i 19841) graphite com pared to iHanneJ il988l) amorphous carbon and to 
iDraine & Legh984 astronomical silicate. The coefficients were computed 
using Mie scattering routines using an MRN size distribution with mini- 
mum and maximum grain radii of 0.005 and 0.05 fj,m, identical to that used 
in all models presented in this paper. 



of 0.1 - 2 • 10"^ Mq, again ruling out a significant variation in 
dust mass between day 775 and day 6067. Our own dust mass es- 
timates in Table [T]imply that for the graphite models the total dust 
mass hardly increased between days 615 and 775. Finally, Bouchet 
et al. (2004) detected mid-IR silicate emission from the circumstel- 
lar ring that is believed to have been produced by mass loss by the 
progenitor star of SN 1987 A. This ring is responsible for most of 
the present-day mid-IR emission, but Bouchet et al. (2004) showed 
that this contribution started after day 4200 only, and therefore does 
not affect the present findings. 

We also made estimates for the circumstellar extinction in the 
envelope of SN 1987 A at late epochs in order to constrain our 
m odels. For this w e used the methods described in greater detail 
by ISugerman etal] (l2006h for the case of SN 2003gd. Using the 
photometric da ta of Ham uv & Suntzefl] ([l990') for days 1-813 and 
ICaldwellet ail (1993) thereafter, we estimated the circumstellar R- 
band extinction of SN 1987 A at days 615 and 775 by comparing 
the observed level of the light curve at those epochs to that expected 
for a light curve powered purely by the decay of ^®Co that had been 
normalised to the photometric data obtained earlier than day 400. 
This yielded Ar's of 1.7 and 2.4 magnitudes at days 615 and 775, 
respectively. A somewhat more conservative estimate of the cir- 
cumstellar extinction can be obtained if we assume that the intrinsic 
light curve has contributions from other sources of decay, such as 
^^Co, and also allow for the decreasing opacity of the envelope to 
the gamma rays that power the light curve. This yielded yl^ 's of 1 .2 
and 1.6 magnitudes at days 615 and 775, respectively. Correspond- 
ing estimates for the R-band circumstellar extinction produced by 
the models were obtained by comparing the R-band flux levels ob- 
tained with and without dust. Inspection of Table 1 (final row for 
each model) shows that all of our day 615 models satisfy the Ar = 
1.2-1.7 observational constraint. For the case of our day 775 mod- 
els the predicted Ar range is 1.5-2.4, which can be compared to 
our estimated observational constraint of Ar = 1.6-2.4. 

For Clumpy I amorphous carbon or graphite models, the de- 
rived day 615 and day 775 dust masses are the same, while the 
inner radius for the day 775 models is 5.7/5.0 times larger than for 
day 615. Yet the day 775 models yield higher Ar's than the day 



Figure 13. Alternative Clumpy II model fits to the day 615 optical-IR spec- 
trophotometry published by Wooden et al. (1993) -black line- . All mod- 
els used an MRN size distribution truncated to max imum grain radii of 
0.05mm. The red dotted line con'esponds to a graphite j Praine &Lej 19841) 
model. The dashed green line corresponds to a model with a composition of 
50%:50% astronomical silicates : graphite by mass (Draine & Lee 1984]), 
the blue dash-dot line corresponds to 25%:75% astronomical silicates : 
graphite, and the c yan dotted line coiTe sponds to 15%:85% astronomical 
silicates : graphite toraine&Le3ll984 . The mixed composition models 
shown all correspond to a dust mass of 1.3-10~^ Mq, twice as large as for 
the Clumpy II day 615 graphite model that is plotted. 



615 models. The equivalent optical depth increase from day 615 to 
775 is due to the fact that for the volume filling factor to decrease 
from 5 to 3% whilst keeping the same dust mass , the dust number 
density in each individual clump must increase, hence making it 
harder for radiation produced in the internal regions of the clumps 
to escape. The decrease of the clump covering factor from day 615 
to 775 is physically justified under the assumption that no clumps 
and no extra dust condensed between day 615 to day 775. If this 
were true then the expansion of the inner radius from 5 x lO^'^cm to 
5.7 X lO^^cm would imply a decrease of the volume covering factor 
from / = 0.05 to / = 0.034. However, best fit of the SED implies / = 
0.03, causing therefore the increase of Ar from 1.4 to 2.0 reported 
in Table 1 . 

The behaviour of Clumpy II models is somewhat different; the 
volume filling factor here has to increase from 0.1 to 0.3 in order 
to fit the decrease of the SED in the optical from day 615 to day 
775. As discussed in Section 3.2, once the clumps become optically 
thick, the only factor affecting the optical SED is /, with the main 
effect of varying the dust mass being on the temperature of the 
grains and therefore on the IR emission. In this model the total dust 
mass also increases between the two epochs, implying that as new 
grains form new clumps also form. An alternative interpretation is 
that clumps form early on by Rayleigh-Taylor instabilities in the 
explosion, and remain imprinted in the ejecta. The increase in dust 
mass, might therefore be due to additional accretion. 

Our dust mass estimates are dependent on the assumed dust 
chemistry. As shown in Figure [121 amorphous carbon grains with 
Hanner (1998) optical constants have higher absorption coefficients 
per unit mass in the optical an d near-IR than either astronomical 
graphite or silicate grains with iDraine & Le3 ( Il984t) optical con- 
stants. The amorphous carbon models therefore require the least 
mass to fit the observed SEDs of SN 1987A. Graphite models re- 
quire between 2 and 3 times as much mass to provide a fit. Silicate- 
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Figure 14. Alternative Clumpy II models (as in Figure [13) superimposed 
to the 7-12 mm region of the observed day 615 spe ctrum from Figure 1 
of W93. The red dotted line coiresponds to a graphite jPraine & Leel 19841) 
model. The dashed green line coiresponds to a model with a composition o f 
50%:50% astronomical silicates : graphite by mass jPraine &Ledll984) . 
the blue dash-dot line corresponds to 25%:75% astronomical silicates : 
graphite, and the c yan dotted line corr esponds to 15%:85% astronomical 



graphite, and the c yan dotted Ime corr t 
silicates : graphite jPraine & Leell984l) 



only models would require even larger dust masses, due to their 
even lower optical and near-IR absorption coefficients; however 
pure silicate models are ruled out for SN 1987A by the lack of 
detectable 10- and 18-/im silicate emission or absorption features 
in any of the mid-IR spectra of SN 1987A (W93; Roche et al. 1993; 
Bouchet & Danziger 1993). 

After day 530, the emission-line strengths of [Si I] \.()5^,m, 
[Mg I] 4571 A, and [O I] 6300A decreased significantly ("Lucy et al.' 
[l99 1; Bouchet 8l Danziger 1993), consistent with the formation of 
some fraction of silicate dust. We therefore tried a number of mix- 
tures of astronomical silicates and graphite. Figure [T3l shows fits 
to the day 615 W93 data obtained with our Clumpy II model us- 
ing 50% : 50% (green dash-dot line), 25% : 75% (blue dash-dot 
line) and 15% : 85% (cyan dash-dot-dot line) astronomical sili- 
cate : graphite mixtures, all with masses of 1.3-1O~'^M0 (red dot- 
ted line), i.e. twice as large as for the Clumpy II graphite model. 
The 50% : 50% silicate:graphite and 25% : 75% mixtures produce 
too strong 10- and 18-/im silicate absorption features compared to 
the observational limits. However, the 15% : 85% silicate: graphite 
model produces a smooth enough mid-IR spectrum to be consistent 
with the observations. 

Due to the high optical absorption coefficients of amorphous 
carbon grains, the Clumpy II day 615 amorphous carbon mass esti- 
mate of 2.2-10"^ M0 (Table[TJ is likely to represent a lower limit 
to the total mass of dust that condensed in the ejecta of SN 1987A, 
while the mixed silicate:graphite mass of 1.3-10"^ M0, discussed 
in the previous paragraph, probably represents close to an upper 
limit, h'on grains have optical absorption coefficients in the opti- 
cal that are intermediate between those of astronomical graphite 
and silicate, while lacking the strong infrared features that elimi- 
nate silicate-only models. However, direct observational evidence 
for iron grains in circumstellar or interstellar environments is still 
lacking. 

Figure[T5]shows a comparison between th e day 678 optical-IR 
spectrum of SN 2003gd jSugerman et al.l2006h and the day 615 and 
day 775 spectra of SN 1987 A (W93), after scaling to a common dis- 



Figure 15. The observed SED for day 678 of SN 2003gd (asterisks; Sug- 
erman et al. 2006, and references therein) compared to those of SN 1987A 
on day 615 (dotted line) and day 775 (dashed line). The optical to mid- 
IR spectrophotometry of SN 1987A is from W93 and the 50- and 100 -/jm 
photometric points at day 793 (triangles) are from lHarvev et alj Jl989l) . 



tance of 5 1 kpc. It shows similarities between their SEDs in the op- 
tical and I6-30-/im regions, but difference s in the 3-8-/xm region . 
The day 678 clumpy model for SN 2003gd jSugerman et al J2006h , 
required a dust mass that is 13 times higher than the maximum 
mass allowed for mixed graphite:silicate models of SN I987A and 
27 times higher than our day 775 Clumpy II pure graphite mass 
estimate for SN 1987A. The main reason for this difference is the 
very different opacities of the dust mixtures required to fit the SEDs 
of the two SNe; in the case of SN 2003gd, a dust mixture com- 
posed of 15% amorphous carbon and 85% astronomical silicates 
provided the best fit ( Sugerman et al. 2006), while the maximum 
silicate mass fraction allowed for SN I987A is only 15%, due to 
its featureless mid-IR spectrum at late epochs. Astronomical sil- 
icates are much more transparent to optical and near-IR radiation 
than are amorphous carbon or graphite grains, and therefore require 
larger dust masses to match a given SED. Furthermore, the much 
higher optical albedos of silicate grains allows much larger masses 
to be inside clumps, for similar levels of IR emission, though the 
silicate 10-/xm feature should be present in emission or absorption 
(mainly the latter). We note that models containing only carbon 
grains did not provide a good fit to the day 499 optical-IR spec- 
trum of SN 2003gd, producing too much emission in the 2~10-/im 
region. 

The model fits imply that the proportions of silicates and car- 
bon grains were reversed in SN 1987 A and SN 2003gd. The mass of 
the immedi ate progenitor of SN 2003gd is es t imated to have bee n 
8I2 M o i lVan Dvk. Li, & Filippenkol |2003| ; ISmartt et al.l l2004h , 
while the mass of Sk-69 202, the prog enitor of SN 1987 A, is esti- 
mated to have been ~20 M© jAmett et al. 1989). Nucleosynthetic 
calculations show that although large masses of oxygen should be 
produced by high mass supemovae, below 20 M the yield of 
oxygen from the O-rich z one can drop below the yield of car- 
bon from the C-rich zone jAmett et alj|l989l ; IWoosIev & Weaver! 
1995), although the results are highly sensit ive to details of some 
of th e adopted reaction rates |Arnettl I996h. Irrespective of this, 
Clavt on. Liu & Dalgaiiiol h999l) and iDeneault. Clayton & Hegeil 
( 2006)) have shown that due to the destruction by ionizing particles 
of CO, and the consequent disruption of the normal CO sink for 
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carbon atoms in oxygen-rich regions, large quantities of carbon- 
rich dust particles can form easily in SN ejecta zones that have 
C/0<1. What remains to be elucidated is what causes either sili- 
cate or carbon grains to dominate the dust produced by particular 
supernovae. 



5 DISCUSSION 

We have studied the effects of a number of clumping parameters on 
the SEDs emerging from two geometrical configurations, consist- 
ing of a spherical shell of clumpy medium illuminated by a diffuse 
radiation source. Clumpy I models assume full mixing between the 
dust and the radiating medium, while Clumpy 11 models assume 
that dust is only present in the clumps and the illuminating radia- 
tion comes only from the ICM. 

The main conclusions that can can be drawn from our param- 
eter investigation can be summarised as follows: 

(i) varying the volume filling factor, / has opposite effects in 
Clumpy I and Clumpy II models. For a constant dust mass, Md and 
clump to ICM ratio, a, an increase in / in Clumpy I models causes 
a decrease in the effective optical depth, with less optical radia- 
tion being reprocessed at IR wavelengths. In Clumpy II models / 
directly controls the amount of optical radiation being intercepted 
by the dust, since the illumination source is always external of the 
clumps. This is the same behaviour shown by centrally illuminated 
shell models. 

(ii) in the range explored, variations of a in Clumpy I models 
(constant A/d and constant /) have only a minor impact on the 
emerging SEDs. 

(iii) for models where the emitting source is intermixed with 
the dusty medium, dust grains cannot be hidden in optically thick 
clumps. For these models a degeneracy exists between i\/d and / 
which may not be resolved with IR SED observations alone. 

(iv) for Clumpy II models with optically thick clumps, an in- 
crease in Md results in cooler grains which will therefore radiate at 
longer IR wavelengths. 

We used both smoothly varying and clumped dust density dis- 
tributions to obtain new estimates for the mass of dust condensed 
by the Typ e II SN 1987A by fitti ng the 0.3-30^im spectrophotomet- 
ric data of I Wooden et alj ( Il993h at two epochs (day 615 and day 
775). From our modelling we conclude that between 2x 10~* and 
1.3x10^'' Mq of mainly carbon-based grains had formed in the 
ejecta by day 615 (4-7 x 10^^ Mq) for pure graphite models) with 
the derived dust masses being found to be independent of whether 
the dust was clumped or not. This result differs from the conclu- 
sions of W93, who derived unclumped graphite dust masses of 4- 
6xlO~^ Mq on days 615 and 775, versus graphite masses from 
their analytic clumped dust modelling of 3-5 x 10"^ M0 . The latter 
values are close to our own clumped graphite dust mass estimates 
from Monte Carlo modelling. For a 20 Mq SN progenitor, between 
0.3-0.4 Mq of refractory elements are expected to form (Woosley 

6 Weaver 1995), so the estimate of 2x10"" to 1.3x10"^ Mq 
of new dust formed in the ejecta of SN 1987A corresponds to a 
condensation efficiency of only 5x 10~" to 4x 10~^, significantly 
lower than the dust condensation efficiency of ~0. 12 estimated for 
SN 2003gd by Sugerman et al. (2006). The cause of such very 
different dust yields amongst Type II core-collapse supernovae re- 
mains to be explained. 
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